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Abstract 

The scenario of isospin-violating dark matter (IVDM) with destructive interference between DM-proton and DM- 
neutron scatterings provides a potential possibility to reconcile the experimental results of DAMA, CoGeNT and 
XENON. We explore the constraints on the IVDM from other direct detection experiments such as CRESST and 
SIMPLE, etc. and from the indirect DM searches such as the antiproton flux measured by BESS-Polar II. The results 
show that the relevant couplings in IVDM scenario are severely constrained. 
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Some of the recent dark matter (DM) direct detec- 
tion experiments such as DAMA G] |2l [3, CoGeNT 
El [51 and CRESST-II |6l have reported events which 
cannot be explained by conventional backgrounds. The 
excesses, if interpreted in terms of DM particle elastic 
scattering off target nuclei, may imply light DM parti- 
cles with mass around 8-10 GeV and scattering cross 
section around 10 cm^. Other experiments such as 
CDMS-II cm, XENONlO/100 lElIIOl, and SIMPLE 
flTI etc., have reported null results in the same DM 
mass range. 

A commonly adopted assumption on interpreting the 
DM direct detection data is that in spin-independent 
scatterings the DM particle couplings to proton (fp) and 
to neutron (/„) are nearly the same, i.e. /„ x fp, which 
makes it straight forward to extract the DM-nucleon 
scattering cross sections. It is a good approximation 
for neutralino DM and DM models with Higgs portal, 
e.g, the scalar DM in left-right models fTIl [131 [T5l 
and 4th generation Majorana neutrino DM lfT6l . How- 
ever, in generic cases, the interactions may be isospin- 
violating fill [181 [191 EOl Ell |22l|23]|. In this scenario, 
the DM particle couples to proton and neutron with dif- 
ferent strengths, possible destructive interference be- 



tween the two coupUngs can weaken the bounds of 
XENON 10/ 100 and move the signal regions of DAMA 
and CoGeNT to be closer to each other (TP, '22]. In 
order to reconcile the data of DAMA, CoGeNT and 
XENONlO/100, a large destructive interference corre- 
sponding to fn/fp ~ -0.7 is required f^Tl. 

Possible constraints on IVDM from the cosmic neu- 
trinos and gamma ray on IVDM have been discussed 
previously in Refs. lH |25l EH . Recently the BESS- 
Polar II experiment has measured the antiproton flux 
in the energy range from 0.2 GeV to 3.5 GeV lIZTl 
which have higher precision compared with that from 
PAMELA 1 28 1 at low energies. In this talk, we discuss 
on the direct and indirect constraints on IVDM with fo- 
cus on the cosmic -ray antiproton constraints. The de- 
tails of our analysis can be found in Ref. [29 |. 

For a DM particle x elastically scattering off a target 
nucleus, the differential scattering cross section can be 
written as 



dcr 
dER 



-o"o , 



(1) 



where F(Er) is the form factor of the nucleon and 
Ha - {fn^fnAilim^ + m^) is the DM-nucleus reduced 
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mass. The quantity ctq can be understood as the total 
scattering cross section at the limit of zero-momentum 
transfer which is related to /p(„) through 



cro = ^[z/^ + (A-Z)/„]' 



(2) 



where Z is the atomic number and A is the atomic mass 
number. Under the assumption that the scattering is 
isospin conserving (IC), i.e., /„ « fp, the total cross 
section ctq is independent of Z and only proportional to 
A^. One can define a cross section cr'^ which is the 
value of (Tp extracted from ctq under the assumption of 
IC interaction as 



CTQ . 



(3) 



In the generic case where /„ + fp, the true value of cr^ 
will differ from cr'^ by a factor F{f„/fp) which depends 
on the ratio /„/ fp and the target material 



o-p = F(f„/fp)cr 



(4) 



If the target material consists of kind of relevant nu- 
clei with atomic numbers {a = 1 , . . . , A^) and frac- 
tional number abundances Ka, and for each nucleus Z„ 
there exists M type of isotopes found in nature with 
atomic mass number Aq,,- and fractional number abun- 
dance rjai (/ = 1, . . . , M), the expression of F{ f„/fp) can 
be explicitly written as 



F(f„/fp) = 



(5) 



where jua^i is the reduced mass for the DM and the nu- 
cleus with atomic mass number A^i- For a given target 
material T, there is a particular value of f„/fp which 
corresponds to the maximal possible value of F{f,J fp) 



(6) 



The value of varies with target material. In Tab. [T] 
we list the values of for some typical material utilized 
by the current or future experiments. 

If the values of the target material used by two 
experiments are very close to each other, the tension be- 
tween the two experimental results, if exists, is less af- 
fected by the effect of isospin violation. From Tab. [T] 
one finds that ^Na ~ fcciF, = -0.92, ^xe ~ ^Csi ~ -0.7 
and ^si ~ ^Ca(W)04 - -1-0. Thus the tension between 



DAMA signal from Na recoil and the upper bound from 
SIMPLE is unlikely to be alleviated by isospin viola- 
tion, which can be clearly seen in Fig. [T] Similarly, if 
there exists contradictions between XENON and KIMS, 
CoGeNT and the Ar based experiments such as Dark- 
Side, it can hardly be explained by isospin violating 
scattering. The SIMPLE result is also useful in compar- 
ing with the CRESST-II which utilizes Ca(W)04 which 
has ^ca(W)04 = -1.0. Obviously, for the experiments 
use the same target material, the possible tension be- 
tween them cannot be relaxed by isospin violation, such 
as the tension between CoGeNT and CDMS-II, as both 
use germanium as target nucleus. 

In Fig. [T] the allowed regions by the current experi- 
ments are shown in the {(jp, m^) plane for f„l fp - -0.70 
. It can be seen that the overlapping region between 
GoGeNT and DAMA may still be consistent with the 
exclusion curve from the XENON 100 2011 data [fOl. 
However, If one considers the recently updated upper 
bounds from XENON 100 |30|, the main bulk of the 
overlapping region is excluded for both the GoGeNT 
results with and without surface event rejection correc- 
tions, which challenges the IVDM as a scenario to rec- 
oncile the results of DAMA, CoGeNT and XENON. 
The overlapping region between DAMA and CoGeNT 
seems also to be excluded by the results of SIMPLE 
1,11 J and CDMS-II independently |i7j|8J. Note however 
that there still exists controversies regarding the detec- 
tor stability of SIMPLE experiments 1311 ES, the re- 
coil energy calibration of CDMS experiment |(33l| and 
the extrapolation of the measured scintillation efficiency 
to lower recoil energy in the previous XENONIOO data 
analysis 1,34.. 35 J . 

We assume that the DM particles interact with the 
SM light quarks through some heavy mediator particles 
much heavier than the DM particle such that both the 
scattering and the annihilation processes can be effec- 
tively described by a set of high dimensional contact 
operators 



X = ^ aiqOiq . 



(7) 



If the DM particles are Dirac fermions, the relevant op- 
erators arising from scalar or pseudoscalar interactions 
are given by 



0\q = xxm^Oiq^xj^xm, 
03q = xxqr'q^OAq^xT^xqT^q- 



(8) 



The operators from vector or axial-vector type interac- 
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Xe 


Ge 


Si 


Na(I) 


Ca(W)04 


C2CIF5 


Csl 


Ar 


-0.70 


-0.79 


-1.0 


-0.92(-0.73) 


-1.0(-0.69) 


-0.92 


-0.71 


-0.82 



Table 1 : Values of ^7 for different target material. For Nal, the two values -0.92 and -0.73 correspond to the scattering olf Na and Nal respectively. 
Similarly, for CaW04, the two values -1.0 and -0.69 corresponds to the scattering without and with tungsten nuclei respectively. 
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Figure 1: Favored regions and limits in the {crp,my) plane from var- 
ious experiments for /;,//;, = -0.70 such as DAMA |2|, GoGeNT 
unmodulated data |5|, CoGeNT unmodulated data with surface 
event rejection factors taken from Kelso, etal |36|, CRESST-II [6J, 
XENONlO/100 19 301, CDMS EUS and SIMPLE jlT). 



tions are 

Oiq = xy''xq7iiy^q, Osq = xyy^xqyi^r^q, (9) 

and the ones from the tensor interactions are 

Ogq ^ xcr^^'xqo-fivq, Owq ^ x^^^y^xq^-^vq- (lO) 

If the DM particles are Majorana particles the vector 
and tensor operators are vanishing identically. Among 
these operators only Oiq and 0^^ contribute to spin- 
independent scattering cross sections at low velocities. 
The scattering cross sections induced by the operators 
02g and 0(,q are velocity suppressed. The operators 
Ojg and Osq contribute only to spin-dependent scatter- 
ing cross section, and the nucleus matrix elements for 
the operators O^q, O^q, Ogq and Ou)q are either vanishing 
or negligible. The DM annihilation into quarks through 
Oiq is a p-wave process, which is velocity suppressed. 
It does not contribute to the cosmic antiproton flux, but 
still contributes to the DM relic density as p-wave pro- 
cesses is non-negligible at freeze out. 



Similarly, for DM being a complex scalar <p, possible 
operators are 

On - 2m^(4)*cf>)qq, On - 2m^((f>*(f))qy^q, 

On = (cl,%4,)qy'q,0u^(^%cl,)qy''^in) 

Among those only 0\\ and On contribute to the spin- 
independent scatterings. The DM annihilations through 
operator On are p-wave processes. 

In summary, we only consider four operators 



Olq, Osq, O 



llq: 



andO 



13q, 



which are relevant to IVDM. 

The DM couplings to nucleons fp „ can be expressed 
in terms of the DM couplings to quarks follows 



(12) 



For the Dirac DM with scalar interaction aujxxqq, one 
has Bf*"' = Z^^"' mp(„)/mq for q - u,d,s and B'^^"'' - 



\q JTq 



fP(n) 



is the 



(2 /27) fj^ mp(„)/mq for q - c,b,t, where fj^^ 
DM coupling to light quarks obtained from the cr-term 
{N\mq-qq\N) = /^M^, and f^^^ = 1 - 1;,=,,,,^,., /^f ■ In 
order to maximize the isospin violating eff'ect, the co- 
efficients B^fi^ 1^ J, must be strongly suppressed. As- 
suming that the DM-nucleon couplings are dominated 
by the DM couplings to the first generation quarks, the 
ratio /„//p is given by 



(13) 



The value of f„/fp - -0.7 can be translated into 
a\dla\u - -0.93 at quark level. This value is the same 
for complex scalar DM. For operator Osq one simply 
has B^f' = 2(1) and B^f = 1(2), and fif " = for 
q = c,s, t, b. 



fn_ _ flSu + 2a5j 

fp 2fl5„ + a^d 



(14) 



Thus for /„/ fp - -0.7, one finds a^d/asu - -0.89. 

Annihilation or decay of light DM particles in the 
galactic halo can contribute to exotic primary sources 
of the low energy cosmic ray antiprotons, which can 
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be probed or constrained by the current satellite- 
and balloon-borne experiments such as PAMELA and 
BESS-polar II, etc.. The predicted antiproton flux from 
DM annihilation depends on models of the cosmic-ray 
transportation, the distribution of Galactic gas, radi- 
ation field and magnetic field, etc.. It also depends 
on the particle and nuclear interaction cross sections. 
In this work, we use the numerical code GALPROP 
Il371 ESI [m EOl m] which utilizes realistic astronom- 
ical information on the distribution of interstellar gas 
and other data as input and consider various kinds of 
data including primary and secondary nuclei, electrons 
and positrons, y-rays, synchrotron and radiation etc. in 
a self-consistent way. In the GALPROP approach, we 
consider several diffusion models (parameter configura- 
tions). The different results between the models can be 
regarded as an estimate of theoretical uncertainties. 

In the diffusion models of cosmic ray propagation, 
the Galactic halo where diffusion occurs is parame- 
terized by a cylinder with half height Z/, and radius 

= 20 - 30 kpc. The densities of cosmic ray parti- 
cles are vanishing at the boundary of the halo. The pro- 
cesses of energy losses, reacceleration and annihilation 
take place in the Galactic disc. The source terms for 
the secondary cosmic rays are also confined within the 
disc. The diffusion equation for the cosmic ray particle 
is given by 



i - ^^^^ 

_d_ 



- — l/r- —il/ + q{r,p), 
Tf Tr 



(15) 



where i//{r, p,t) is the number density per unit of to- 
tal particle momentum which is related to the phase 
space density /(r, p, t) as i/'(r, p, t) - Anp^fir, p, t). For 
steady-state diffusion, it is assumed that di/z/dt - 0. The 
spatial diffusion coefficient D^^ is parameterized as 



-./3Dq 



(16) 



where p = p/{Ze) is the rigidity of cosmic ray parti- 
cle and 6 is the power spectral index which may take 
different values or ^2 when p is below or above the 
reference rigidity po- is a normalization constant, 
and p - v/c is the velocity of the cosmic ray particle. 
The convection term is related to the drift of antiproton 
from the Galactic disc due to the Galactic wind. The 
direction of the wind is usually assumed to be along the 
z-direction which is perpendicular to the galactic disc 



and is a constant Vc = [20(1) - l]Vc. The diffusion in 
momentum space is described by the reacceleration pa- 
rameter Dpp which is related to the Alfven speed Va of 
the disturbances in the hydrodynamical plasma as 



2„2 



(17) 



3D,,5(4 - 6^)(4 - 6)w' 

where w stands for the level of turbulence which can be 
taken as w = 1 as Dpp depends only on the combination 
V^/w. p is related to the momentum loss rate, t/ and 
Tr are the time scale for fragmentation and radioactive 
decay respectively. 

Thus the whole diffusion process depends on a num- 
ber of parameters: R, Zi,, po. Do, 61/62, Vc , Va, Ps and 
yi/ji- The recently updated analysis based on Markov 
Chain Monte Carlo fits to the astrophysical data show 
that Zi, should be around 4 - 7 kpc ||42l 03]. In or- 
der to obtain a conservative upper bounds we choose 
Z/i a! 4 kpc in our analysis. 

We consider several typical propagation models in 
GALPROP, and focus on the models with the secondary 
antiproton background below the current data, which 
leaves room for DM contribution and results in conser- 
vative upper bounds. The first one is the plain diffu- 
sion model (referred to as "Plain") in which there is no 
reacceleration term f4T\. The second one is the conven- 
tional model (referred to as "Conventional") with reac- 
celeration included 1391 1411 . The last one (referred to 
as "Global-Fit") is the model from a global fit to the 
relevant astrophysical observables using Markov-Chain 
Monte Carlo method ll43l . The main parameters of the 
three models are listed in Tab. |2] 

The primary source term from the DM annihilation 
has the form 



q(r) = 77n(r)^<crvrei) 



dN_ 
dp ' 



(18) 



where n(r) - p(r)/m^, r] - 1/2(1/4) if the DM particle 
is (not) its own antiparticle and dN/dp is the injection 
spectrum per DM annihilation. For the DM profile we 
took the isothermal profile 1441 



P(r) = po 



r^+Rl 



(19) 



where = 8.5 kpc is the distance of Solar system from 
the galactic center, Rs - 2.8 kpc and the local den- 
sity is taken to be po = 0.3GeVcm"-'. The choice 
of isothermal profile and the local density is again to 
achieve a conservative estimate of the antiproton flux, if 
the NFW profile is used, the predicted flux will be en- 
hanced roughly by at most 70%, thus more severe con- 
straints are expected. 
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model 


R(kpc) 


Z/,(kpc) 


Do 


Po 


5lM2 


V«(km/s) 


Ps 


7i/72 


Plain 


30 


4.0 


2.2 


3 


0/0.60 





40 


2.30/2.15 


Conv. 


20 


4.0 


5.75 


4 


0.34/0.34 


36 


9 


1.82/2.36 


Fit 


20 


3.9 


6.59 


4 


0.3/0.3 


39.2 


10 


1.91/2.40 



Table 2: Propagation parameters in the "Plain" 1411 , "Conventional" I39l|41| and "Global-Fit" [43] models used in the GALPROP code. Do is in 
units of lO^^cm^ ■ s"', the break rigidities po and p, are in units of GV. 



In the case where the DM particles are thermal relics, 
possible large couplings to light quarks may under pre- 
dict the DM relic abundance in comparison with the ob- 
served value = 0.113 ± 0.004 |45|. Thus the relic 
density can also impose upper bounds on the relevant 
couplings. Whereas the p-wave annihilation processes 
give no contribution to the antiproton flux, in the cal- 
culation of DM relic density, the p-wave annihilation 
is nonnegligible as the typical relative velocity of DM 
particles is Vrei/c ~ (9(0.3) at freeze out. The annihi- 
lation cross section times the relative velocity can be 
expressed as crviei = a + bv^^^, where a and b are co- 
efficients corresponding to the .s-wave and p-wave con- 
tributions. For light DM around 10 GeV, - 61.75, 
we find Xf » 22. The relic abundance is approximately 
given by 

1.07 X 10^ GeV- ^^^^ 
Mpi yjgl a + 3bxf 

We perform analysis to obtain upper bounds on the 
DM isospin violating couplings to light quarks. For the 
sake of simplicity, we consider one operator at a time 
and ignore interference between the operators. 

For the operators contribute to i-wave annihilation 
which leads to velocity-independent annihilation cross 
sections, the relevant DM couplings to quarks are found 
to be tightly constrained by both the cosmic antipro- 
ton flux and the thermal relic density. In Fig. |2] the 
constraints on the coeflicients a^g for operator Osq - 
Xy'^X'iJiil are shown in the (asu,asd/a5u) plane. The 
mass of the Dirac DM particle is fixed at - 8 GeV. 
The results show that the DAMA- and CoGeNT-favored 
regions are in tension with both the cosmic antiproton 
flux and the thermal relic density. 

At asd/osii - -0.89 which corresponds to f„/fp - 
-0.70, the DAMA and CoGeNT favored value is 05,, ^ 
1.6x 10"^ GeV"^ corresponding to an annihilation cross 
section of (crviei) ^ 3.4 x 10""^ cm^. As it can be 
seen in Fig. [3] for such a large cross section, the pre- 
dicted antiproton flux is much higher than the current 
BESS-Polar II and PAMELA data and results in a huge 
;^rVd.o.f=1.3 X IOV35 in the "Global-Fit" model and 



DAMA 3a 




Figure 2: Upper bounds on the coefficient asu as a function of asj /asu 
at 95% CL from cosmic antiproton flux and DM relic density. The 
mass of DM particle is fixed at 8 GeV. The favored regions and exclu- 
sion contours from various experiments such as DAMA 12 ], GoGeNT 
(5J, XENON |9ji30J, CDMS tlji^i and SIMPLE lllj are also shown. 

an even larger one in the "Conventional" model. The 
upper bound set by the antiproton data at 95% CL is 
a5„ < 1.7 X 10 GeW^ in the "Global-Fit " model at 
dsd/cisu = -0.89, which is about an order of magni- 
tude lower and corresponds to an annihilation cross sec- 
tion (crvjei) ~ 3.7 X lO^^^cm^. The predicted antiproton 
fluxes are also shown in Fig. |3] In the two propaga- 
tion models, the upper bound from the "Conventional" 
model is slightly stronger than that from the "Global- 
Fit" model. In Fig. |3] we also plot the antiproton back- 
ground of "Plain" model which is already higher than 
the data. Thus the upper bounds from this propaga- 
tion model are expected to be much stronger. Since 
we are interested in conservative upper bounds, we do 
not further investigate the constraints in this model. 
In Fig. [3] the upper bound from the relic density is 
asu < 6.0 X 10"^ GeV"^ at a^d/asu = -0.89, which is 
weaker than that from the antiproton flux but still in ten- 
sion with the DAMA- and CoGeNT-favored value. 

Similar results are found for the operator Onq ll29l . 
Compared with the case of Osg, the constraints on the 
coefficients of Ong are weaker, which is due to the fact 
that for the hadronic matrix element of scalar operator 
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Figure 3: Predictions of cosmic antiproton spectra from DM an- 
nihilation induced by operator Osq in tlie "Global-Fit" propagation 
model. Two cases are considered: (A) For as,, = 1.6 X 10"'' GeV"^ 
which is favored by the DAM A and CoGeNT experiments. (B) For 
i^Su = 1.7 X 10"'' GeV"- which is the maximal value allowed by the 
cosmic antiproton data at 95% CL. The ratio asj/asi, is fixed at -0.89 
corresponding to /„//,, = -0.70 and the mass of DM particle is fixed 
at 8 GeV. The data of BESS-Polar II 1 27 1 and PAMELA 1 28 1 are also 
shown; right ) Same as left), but for the "Conventional" model. 



qq the B,v, factors are larger than that for vector operator 
qy'^q, which allows smaller a,-^ for the same value of fp^„ 
and results in smaller annihilation cross sections. The 
constraint from the thermal relic density is in tension 
with the DAMA- and CoGeNT-favored values. The op- 

erators Oiq - XXII ^"d Ong = (0* 5^ 4')qyq contribute 
to p-wave annihilation with cross section proportional 
to V 



reV 



Thus they do contribute very little to the cos- 
mic antiproton flux. However, their contributions to the 
thermal relic density cannot be neglected, as at freeze 
out the relative velocity is finite. For the operator Oi^q 
one can see some tension between bounds set by the 
relic density and regions favored by DAMA and Co- 
GeNT data. The constraint is not as stringent as that 
from the latest XENONIOO data. For the operator Oig, 
the constraints from relic density is rather weak. The 
difference is again due to the different Bjg factors for 
these two type of interactions. 

It is straight forward to extend the discussions to Ma- 
jorana fermions and real scalars. For the particles being 
its own antiparticles the primary source of the anitpro- 



ton will be enhanced by a factor of 2 in Eq. ( 18 i, which 
may lead to more stringent constraints. 

In summary, we have investigated the allowed val- 
ues of DM-nucleon couplings in the scenario of IVDM 
for various target nuclei used in DM direct detections. 
We find that the recently updated XENONIOO result 
excludes the main bulk of the overlapping signal re- 
gion between DAMA and CoGeNT. We have shown that 



whereas the effect of isospin violating scattering can re- 
lax the tensions between the data of DAMA, CoGeNT 
and XENON, the possible disagreement between some 
group of experiments such as that between DAMA and 
SIMPLE are not likely to be affected for any value of 
fnlfp- We have investigated the conservative constraints 
on the couplings between the IVDM and the SM light 
quarks from the recent cosmic ray antiproton data and 
that from the thermal relic density. Among the four op- 
erators relevant to IVDM Oiq, O^q, Owq, O^q, the oper- 
ators Osq and 0\ \q are found to be tightly constrained by 
the antiproton data and O^q is constrained by the relic 
density. Only the operator 0\q can survive both the con- 
straints while contribute to large enough isospin violat- 
ing interaction required by the current data of DAMA, 
CoGeNT and XENON. 
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